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' The dissipative dynamics of a vortex line in a superfluid is investigated within the frame of a non- 

Markovian quantal Brownian motion model. Our starting point is a recently proposed interaction 
Hamiltonian between the vortex and the superfluid quasiparticle excitations, which is generalized 
to incorporate the effect of scattering from fermion impurities ( 3 He atoms). Thus, a non-Markovian 
equation of motion for the mean value of the vortex position operator is derived within a weak- 
. coupling approximation. Such an equation is shown to yield, in the Markovian and elastic scattering 

' limits, a 3 He contribution to the longitudinal friction coefficient equivalent to that arising from the 

Rayfield-Reif formula. Simultaneous Markov and elastic scattering limits are found, however, to be 
incompatible, since an unexpected breakdown of the Markovian approximation is detected at low 
i-C ' cyclotron frequencies. Then, a non-Markovian expression for the longitudinal friction coefficient 

Jr? | is derived and computed as a function of temperature and 3 He concentration. Such calculations 

S show that cyclotron frequencies within the range 0.01—0.03 ps _1 yield a very good agreement to the 
longitudinal friction figures computed from the Iordanskii and Rayfield-Reif formulas for pure 4 He, 
! 1 1 ' up to temperatures near 1 K. A similar performance is found for nonvanishing 3 He concentrations, 

where the comparison is also shown to be very favorable with respect to the available experimental 

■ data. Memory effects are shown to be weak and increasing with temperature and concentration. 

PACS numbers: 67.40.Vs, 67.40.-w, 67.60.-g, 05.40.Jc 

a- 

^ ' I. INTRODUCTION 

a : 

It is widely accepted that the superfluid vortex dynamics at zero temperature is ruled by the Magnus forced 

m„v = p s KZx(v-v s ). (1.1) 

Here we are assuming a straight vortex line parallel to the z axis, moving with a velocity v. v s denotes a uniform 
background superfluid velocity which, in the simplest case, may be assumed to be time-independent and then dropped 

■ from Eq. ijl.lfl , if v is reinterpreted by Galilean invariance as a vortex velocity relative to a background superfluid at 
y—i \ rest. p s denotes the superfluid mass density and k the quantized circulation of the vortex velocity field (e.g., n — h/m^ 
t-H ■ for one quantum of counterclockwise circulation, being 7714 the mass of a 4 He atom and h the Planck's constant). Then, 

the right-hand side of Eq. (|1.1|) represents the Magnus force per unit length acting on the vortex and, accordingly, 
m v on the left-hand side represents an effective vortex mass per unit length. However, there is no consensus in the 
literature as regards the value of m v . Most of the treatments so far, have neglected m v by assuming that it must 
be equivalent to the hydrodynamic mass of a core of atomic dimensions. 1 Then, from i|l.lfl we get the well-known 
law of zero temperature vortex dynamics, which states that a vortex must move at the velocity which the superfluid 
possesses at the location of the vortex itself^ On the other hand, more recent theories^*^ claim that the vortex mass 
should not be ignored, since it is shown to be logarithmically divergent with the system size, thus exceeding by far the 
core mass. Such a large mass, however, can be shown to be consistent with the above law of vortex dynamics, if the 
O ■ dissipative mechanisms acting at zero temperature are taken into account. In fact, one must consider the effect of the 
vortex coupling to the superfluid which should give rise to dissipation in the form of phonon emission, in close analogy 
^ . to the photon radiation mechanism stemming from an accelerated charge in electrodynamics^^ Another dissipative 
mechanism should arise in ordinary helium from the vortex coupling to the Fermi gas of 3 He impurities. Whatever 
the case, it is remarkable that even for an unbounded system leading to a divergent vortex mass, such dissipative 

■ mechanisms should make the vortex reach the superfluid velocity, in accordance with the fundamental law of zero 
temperature vortex dynamics (see Sec. IIV CJI) . 

At nonvanishing temperatures, in addition to phonon radiation and 3 He scattering, there exists a third dissipative 
mechanism stemming from the vortex scattering of superfluid quasiparticle excitations (phonons and/or rotons). 

Now, let us return to the Eq. (|l.l|l and note that for a background superfluid at rest (v s = 0), the vortex dynamics 
turns out to be identical to the two dimensional one of an electron moving in a uniform magnetic field subjected to 
the Lorentz force. Then, expressing the two dimensional vector v in complex notation as V — v x + iv y , we can rewrite 
(|1.1() in the more compact form 
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v = inv, (1.2) 
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which clearly shows that the vortex will move in a circle at the angular frequency (cyclotron frequency) 

n = ^. (i.3) 

m v 

The above mechanisms of dissipation lead to a complex shift of fi, according to which ljl.2|l becomes 

v = {in cS - Vd )v, (1.4) 

where fi e ff denotes the effective angular frequency into which the unperturbed cyclotron frequency is shifted, and 
Vd > represents a damping frequency that sets the time scale at which V tends to zero. The above equation, 
however, can be written in a more familiar form if we return to the vector notation of Eq. (|1.1|) : 

m v ir = (p s K — D')z x v — Dv, (1-5) 

where 

d 1 = Ps K(i-n cS /n) (i.6) 

D = p s Kv d /n (1.7) 

respectively denote transverse and longitudinal friction coefficients)^ and we have assumed that the normal fluid 
remains at rest. Actually, a vortex in motion may drag the normal fluid in its vicinity, but this effect should be 
negligible below 1 K and we shall restrict our study to such situation^ 

At this point, it is important to notice that the cyclotron motion represented by Eq. 1|1.2|) is also characteristic of 
the helical waves on vortex lines and rings, usually known as Kelvin wavesi&£ In fact, each vortex line element in such 
waves executes motion about the undisturbed line in a circle of radius d and with a frequency w, which approximately 
fulfil 

m v u> 2 d — p s KVi + p s Kuid, (1-8) 

where the amplitude of the deformation d is assumed to be much less than the wavelength A. The above equation 
corresponds to the centripetal component of an expression of the form i|l.lfl . where now v s = — m 9 denotes the local 
self-induced velocity^ which points in a direction opposite to the one of the superfluid velocity field generated by 
the undisturbed vortex line. Therefore, the line velocity in i|l.lfl may be expressed as v = ujdO, where the frequency 
oj will be nonpositive if v points in the same sense as the self-induced velocity. In fact, the solution of the quadratic 
equation 1)1. 8fl yields two frequency branches of opposite sign: 
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whose physical meaning can be easily understood in the limit of long wavelengths (vim v /(p s Kd) << 1, w,; ~ nd/X 2 ). 
That is, the positive fast branch lo+ ~ ft corresponds to the cyclotron motion previously described, while the negative 
slow branch u_ ~ —Vi/d ~ —k/X 2 corresponds to the motion of the vortex element with its local self-induced velocity. 
Then we may see that for a massless vortex line only the slow branch would exist, this being the common assumption 
in the literature of helium vortex waves. As regards experimental studies, only the slow branch has been detected by 
means of a resonant coupling to transverse radio-frequency electric fields acting on vortex lines charged with ions^i 
On the other hand, the thermal excitation of Kelvin waves has been theoretically investigated with rather surprising 
results^ In fact, it was found that the entropy of such waves increases above temperatures about 1.85 K, so that 
the free energy of the vortices is driven negative, with the consequence that superfluidity would be destroyed. This 
phenomenon has been called the "free energy catastrophe" and the authors suggest that it could arise from their 
neglect of the effect of the vortex line on the neighboring phonons and rotons in the system. Now, given that the 
major contribution to the free energy comes from the slow branch, such a "catastrophe" would apparently be shifted 
towards temperatures above the lambda transition if only the excitation of the cyclotron branch were taken into 
account. Actually, this has implicitly been assumed by most of the studies on thermal excitations of vortices through 
phonon and roton scattering, since the authors have only considered straight vortex lines. In particular, the phonon 
scattering excitation of the slow branch was analyzed by Fetter^ and Sonin^ who concluded that it yields only a 
small correction to the friction force calculated for a rectilinear vortex. In fact, the former restriction to considering 
only straight vortex lines in calculations of the friction coefficients, arises naturally if we accept the basic premise 
that it is only the relative motion of an element of line with respect to the normal fluid what matters in such type 
of calculations. That is, any relative motion of a vortex line element should be subjected to the same kind of friction 
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force per unit length, i.e., the same value of the friction coefficients D and D' . Following these considerations, we 
have focused our calculations on the simplest situation of a damped cyclotron motion of a straight vortex line. 

There are, to our knowledge, no experimental results on the transverse friction coefficient below 1.3 K. On the 
other hand, as regards the longitudinal friction coefficient, we must refer to the pioneering experiments performed by 
Rayfield and Reif (R-R) in the early sixties^ In fact, they studied the temperature dependence of the rate of energy 
loss of charge-carrying vortex rings moving through helium II. The radii of such rings are large (> 500 A) compared 
to the distance over which a vortex line is expected to interact appreciably with a quasiparticle. Hence the frictional 
forces on these vortex rings must be the same as those on vortex lines bent into circles. So, R-R were able to measure 
what they called the attenuation coefficient a, which turns out to be simply proportional to the longitudinal friction 
coefficient (a = kD/2). Here it is important to notice that the energy losses in the R-R experiment are consistent 
with a friction owing to the axial displacement of rings, i.e., the main relative motion of each line clement with respect 
to the normal fluid will not correspond to the cyclotron motion. We shall see, however, that in accordance with the 
above basic premise, the longitudinal friction coefficient arising from our theory shows an excellent agreement with 
the one arising from the R-R attenuation coefficient a. We notice also that for a cyclotron motion, the radiation 
damping should be at least comparable to the scattering one for temperatures below 1 K (see Sec. IIV C|l . even though 
we shall focus exclusively on the scattering processes, since phonon emission is supposed to be negligible for the axial 
displacements in the R-R experiment. 

Using kinetic-theory arguments, R-R showed that a comprises three terms, one for each class of quasiparticle 
interacting with the vortex, namely phonons, rotons and 3 He impurities. Each of such contributions was shown to 
be proportional to a corresponding averaged cross section over all momenta, and R-R could determine by fitting to 
their experimental results, that the roton and 3 He cross sections are approximately temperature-independent, with 
respective values 9.5 and 18.3 A. The R-R experiments were performed in the range of temperatures between 0.28 and 
0.7 K, and 3 He concentrations between 1.4xl0~ 7 (ordinary helium) and 2.84xl0 -5 . Then, at the lowest temperatures 
and highest 3 He concentrations, only the 3 He contribution to a is appreciable, allowing its separate study. Analogously, 
in the opposite limit of high temperatures and low 3 He concentrations, only the roton contribution to a survives. 
Unfortunately, only scant information as regards the phonon contribution to a could be derived from such experiments, 
since even though phonon scattering is dominant at the lowest temperatures in pure 4 He, the scattering from 3 He 
impurities becomes the most important contribution in ordinary helium. R-R employed Pitaevskii'gi^ calculation of 
the phonon-scattering cross section to evaluate the phonon contribution to a, but soon after the publication of R-R's 
paper, Iordanskiii& reported an improved theory of the frictional force due to phonons, which seems to be so far the 
most reliable one. Both, Iordanskii's theory and the above kinetic-theory analysis of R-R are based upon an elastic 
scattering assumption, by which the energy of any quasiparticle that collides with the vortex is conserved after the 
collision. This amounts to ignoring any energy the vortex could exchange in such a process, in particular the cyclotron 
energy quantum HQ,, which then should be negligible with respect to the energy of any quasiparticle colliding with the 
vortex. In conclusion, one should expect a cyclotron frequency of finite value, most likely compatible with an elastic 
scattering approximation. Such a hypothesis, has been recently put forward in Ref. ^3 (henceforth to be designated 
as I), where we have studied the friction arising from the scattering of superfluid quasiparticle excitations in the 
form of a translationally invariant interaction potential. Then, the first order expansion in the vortex velocity of 
such a potential was shown to yield vortex transitions between nearest Landau levels, mediated by one-quasiparticle 
transitions. Thus, in the frame of such a model of quantal Brownian motion for the vortex dynamics, the longitudinal 
friction coefficient was computed by making use of weak-coupling and Markov approximations. The result was shown 
to be equivalent, in the limit of elastic scattering, to that arising from the Iordanskii formula and, proposing a simple 
functional form for the scattering amplitude, with a single adjustable parameter whose value was set to get agreement 
to the Iordanskii result for phonons, an excellent agreement with experimental data was found, up to temperatures 
about 1.5 K. Finite values of the cyclotron frequency of order 0.01 ps _1 were also shown to yield practically the same 
results. 

In the present article, we pursue such an investigation in order to analyze the incidence of vortex- 3 He scattering, 
which, as mentioned for ordinary helium, turns out to be the most important contribution to the friction at low 
temperatures. But more importantly, we report an unexpected breakdown of the Markov approximation at low cy- 
clotron frequencies, unnoticed in previous treatments within the elastic scattering limit. Actually, the interaction of 
the vortex with the remaining degrees of freedom of helium, leads to integrodifferential equations of motion for the 
vortex observables, according to which the present vortex motion turns out to be influenced by its whole previous 
history. In the Markov approximation, such a memory is assumed to be negligible and the vortex equations of motion 
are approximated by differential equations like (|1.4|) . 18 We shall show that for low enough cyclotron frequencies, the 
Markov approximation fails and non-Markovian or memory effects must be taken into account. Such effects can be 
of importance in diverse quantum Brownian motion problemsA2*22i2ii2£ and, particularly, in physical situations which 
involve Brownian models of the dynamics of charged particles. For example, a fully non-Markovian reformulation of 
the Abraham-Lorentz theory of radiation reaction in electrodynamics, has been shown to lead to the elimination of 
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"runaway solutions" and causality violations occurring in the original theory* 2 ^ In transport theory, the phenomeno- 
logical Drude-Lorentz result for the ac conductivity has been shown to be affected by important memory effects, 
especially away from resonance 2 ^ and, in the context of two-dimensional magnetotransport, the classical magnetore- 
sistance appears as a consequence of memory effects which are beyond the Boltzmann-Drude approach^S& It may 
be useful to expand on the last problem, since it corresponds just to a classical two-dimensional Brownian motion 
of an electron, subjected to a uniform magnetic field perpendicular to the plane. In fact, the electron is supposed 
to move through a random array of stationary scatterers (background impurities) with short range forces, and there 
are memory effects of two types: (i) the electron may recollide with the same impurity, or (ii) its trajectory may 
repeatedly pass through a space region which is free of impurities. It has recently been shown that backscattering 
processes of the type (ii) are responsible, at low cyclotron frequencies, of additional memory effects leading to unex- 
pected features of the magnetoresistance. 2 ^ Even though there are obviously important differences with the vortex 
problem, it is instructive to compare with this simpler problem, where the source of memory effects at low cyclotron 
frequencies has been fully recognized. 

Microscopic approaches to quantal Brownian motion also show that memory effects are often important when the 
weak-coupling approximation becomes poorer M> This point will be analyzed for our model in Sec. flVBl 

This paper is organized as follows, in the next section, starting from a straightforward generalization to include 
3 He of the Hamiltonian utilized in /, a non-Markovian equation of motion for the vortex dynamics is derived within a 
weak-coupling approximation. Next, we analyze the Markov approximation and the limit of clastic scattering, showing 
that under such approximations, the longitudinal friction coefficient stemming from 3 He scattering, can be shown to 
be equivalent to that arising from the corresponding R-R formula. In Sec. Ill we analyze the breakdown of the 
Markov approximation at low cyclotron frequencies and develop a non-Markovian treatment, from which expressions 
for the longitudinal and transverse friction coefficients are derived. In Sec. IV we focus on the simpler case of a pure 
4 He system. We compare in Sec. IIV Al our results for the longitudinal friction with the Iordanskii (phonon) plus 
the R-R (roton) results, finding a very good agreement within the cyclotron frequency range 0.01—0.03 ps _1 , up to 
temperatures near 1 K. In Sec. IIV Bl we study the frequency ratio f2 e ff/f2, which provides a measure of the memory 
introduced into the vortex dynamics. We explain also the difficulties involved in the calculation of the transverse 
friction coefficient, due to which only its order of magnitude could be estimated. In Sec. IIV CI following the theory 
developed by Arovas and Freire^ we discuss the memory effects related to phonon radiation at zero temperature. 
Finally, Sec. V deals with dilute solutions of 3 He in 4 He, where we compare our results to the available experimental 
data, and extend our study of the memory parameter fl e g/ft in the presence of 3 He. 



II. VORTEX EQUATION OF MOTION, MARKOV APPROXIMATION AND THE LIMIT OF ELASTIC 

SCATTERING 

Our starting point is the following Hamiltonian, which arises as a straightforward generalization of the Hamiltonian 
proposed in /, in order to take into account the presence of 3 He impurities: 

H = H + H lnt , (2.1) 

where 

^ ' k q,cr 

and 

H int = ^ J2 t A ( fc ' #k 5 i + T( - k > (k - q) x z • v. (2.3) 

k,q,cx 

H Q gives the noninteracting part of the Hamiltonian and it comprises three terms, the first of which corresponds to 
the cyclotron motion of the vortex line, the second one to helium II excitations, and the last one to 3 He quasiparticlcs, 
i.e., a^, 6^, and c^ CT respectively denote, a creation operator of right circular quanta, a creation operator of helium 
II quasiparticle excitations of momentum Kk and frequency u>k, and a creation operator of 3 He quasiparticles of 
momentum Tiq, energy e q and spin 1/2 projection a. The interaction Hamiltonian Hi nt arises as a straightforward 
generalization of the form given in /, to include the effect of vortex- 3 He scattering processes. In fact, if we replace 
in Eq. (|2.3(l the vortex velocity operator v as a linear combination of creation and annihilation operators of right 
circular quanta^ it becomes clear that the interaction consists of vortex-quasiparticle scattering events that make the 
vortex to raise or lower one Landau level. Then, in addition to the scattering amplitude A(fc, q) related to the vortex 
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interactions with phonons and rotons discussed in /, now we are including a scattering amplitude q), which takes 
into account vortex- 3 He interactions. 

In previous worka^i^SiSS we derived, by means of a standard reduction-projection procedure and a weak-coupling 
Markov approximation, a generalized master equation for the density operator of the vortex. Our aim was to obtain 
an equation of motion for the mean value of the complex vortex position operator R = x + iy. Now we are interested 
in rederiving such an equation of motion from the more general Hamiltonian Ij2.1|l - H2.3(l . This time we have employed 
a simpler and more direct procedure (see the Appendix), which leads to the following integrodifferential equation of 
motion for v(t) = e^ 24 (R(t)): 

«(*)+/ dTV{r)v{t - t) = 0, (2.4) 

where 







V ^ = ^ o 1 Tl E 4^(k-q) 2 [|A(fc, g )| 2 K-c g )K-n fc ) e '^"^-^ 

+ l\T(k,q)\ 2 (e k -e g )(f q -/^'(WW*-^], (2.5) 

being n k = [exp(?iW/./A;sT) — and f k = {exp[(e/ £ — ^t)/fcsT] + the thermal equilibrium Bose and Fermi 
occupation numbers for the corresponding quasiparticle excitations, respectively. 

The dynamics behind Eq. (|2.4|l can be understood by noting that the present vortex motion is actually influenced 
by its whole previous history, each time being weighted by a memory kernel V(t) such that r — weights the present 
time, while r = t weights the initial condition. Then, if P(r) possesses a microscopic lifetime r m compared to the 
characteristic times that rule the motion of v(t), only the present time t will have a nonnegligible influence upon the 
vortex motion, provided t >> r m . This constitutes the so-called Markov or long time limit approximation^£ii2i22i2ii22 
under which Eq. I|2.4() becomes a differential equation: 

v(t)+uv{t)=0 (2.6) 

where 



/>oo 

v= \ drP(r). (2.7) 
J o 



Recalling that (R{t)) = e int v(t), we may realize that the imaginary part of v yields the shift of the cyclotron frequency 
previously mentioned in Eq. I|1.4|) . i.e., £! c ff = Q — hoav, while the real part, which must be nonnegative, corresponds 
to the damping frequency i>d defined in the same equation. Then, the transverse and longitudinal friction coefficients 
arise from Eqs. (|I.t)|) and 11. 7|) as: 

D' M = (p sK /n)Imv, (2.8) 
D M - ( Ps K/n)Reis, (2.9) 

where the subscript M indicates Markov approximation. The real and imaginary parts of the frequency v, when 
considered as functions of O, obey Kramers-Kronig relations 30 which lead to the following expression for the transverse 
friction coefficient: 

^(fl) = ipf 4,^, (2.10) 

7Ti2 J_ 00 UJ-il 

where P denotes the Cauchy principal part and, 

£>mM = T^jX! ^<z*( k ~ q) 2 [|A(A;, g)| 2 (rz £? - n k )5{w k - uj q - ui) 

^ k,q 

+2|r(fc, q)\ 2 (f q - f k )S(e k /h - ejh - u)]. (2.11) 

The above even function of oj when evaluated at uj — f2 gives the longitudinal friction coefficient, and it is easy to verify 
that Dm{&) > since oj k > uj q =4> n q > n k , and the same for the terms containing the fermion occupation numbers. 
Note that only the scattering events that conserve energy will contribute to the longitudinal friction coefficient (see the 
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arguments of the Dirac deltas in p.llfl for ui = 0). This consequence of the Markov approximation can be physically 
understood in terms of the time-energy uncertainty principle. In fact, in the long time limit only the microscopic 
states with the longest lifetimes are expected to remain with a nonnegligible probability of undergoing a scattering 
transition, and according to the time-energy uncertainty principle, energy should be practically conserved at the end 
of such transitions^ 

We have studied in / the phonon-roton contribution to the longitudinal friction coefficient which arises from the 
first term inside the square brackets in (|2.11|) . We showed that the limit of elastic scattering — > yields an excellent 
agreement with the values derived from experimental data for the roton temperature range, provided the scattering 
amplitude A is set to get agreement with the Iordanskii results for the low-temperature phonon dominated regime. 
Let us now examine the 3 He contribution to the longitudinal friction in Eq. I|2.11|) in the limit f7 — > 0. Denoting such 
a term by D 3 (Q), we have 

d 3 (q) = UmB 8 (n) = -^j£ 4^( k -q) 2 |r(M)l 2 ^( efc - £9 ) 

k,q 

12 



^ Jd 3 k Jd 3 q^ 6(e k ~e q ) (k - q) 2 5(k z - q z ) \T(q, k)\ 2 , (2.12) 



where the continuum limit was explicitly considered in the last expression, A being the area of the system in the x — y 
plane. Here most of the integration can be analytically performed in spherical coordinates, leading to the following 
one-dimensional integral: 

4A 2 m* r^, ilrn ,, l2 df k 



^3(0) = --^- J Q *fc 4 |r(M)r^, (2.13) 

where a Landau-Pomeranchuk dispersion relation e k = h 2 k 2 /2m* was utilized. The above expression can be shown 
to be equivalent to the R-R formula^ provided the scattering amplitude fulfills 

|r(M)l 2 = i^b^( fc M fc )> ( 2 - 14 ) 

where u(k) = hk/m* denotes the group velocity of 3 He quasiparticles and o{k) corresponds to a total momentum- 
transfer cross section for vortex- 3 He scattering, which, in a low 3 He concentration regime^ can be approximated by 
the constant value ao = (18.3 ± 0.7) A. 



III. BREAKDOWN OF THE MARKOV APPROXIMATION. NON-MARKOVIAN TREATMENT 

We have shown in /, that finite values of the cyclotron frequency extracted from recent theories, yield values of 
the longitudinal friction coefficient of the order of that obtained in the elastic limit VI — > 0. We shall henceforth work 
under such an assumption, i.e., Dii(fi) ~ Z?a/(0). Thus it can be shown that Eq. i|2.10|l can be approximated as 
follows: 

2 f°° 

D' M (n) ~ — / dwD M (uj), (3.1) 

where we have also assumed D' M (Vl) >> Dm(Q)- The quasiparticle frequency cutoff in Eq. 1(2. llj) (roughly two 
times the roton frequency) yields the same cutoff to the frequency uj in Eq. I|3.1jl . This shows that the integral in 
(|3.1|) possesses a finite value and thus D' M (n) would diverge as f2 _1 in the limit f2 — > 0. Later we shall see that 
this unphysical result arises from a breakdown of the Markov approximation. In fact, one could expect the effective 
frequency f2 ff to be lower than the cyclotron frequency by the effect of friction, but the existence of a critical cyclotron 
frequency below which the effective frequency becomes negative, seems to be quite unphysical, i.e., one would expect 
effective frequency values bounded as < £l c g = £1 — flD' M (il)/ p s K < ft. We shall see in the following that a non- 
Markovian treatment yields in fact such bounds. To see this, let us return to the integrodifferential equation 12.4(1 
and take its Laplace transform according to the definition v(z) — J Q exp(izt)v(t)dt (Imz > 0): 

V(z)v(z) = izv(z) + u(0), (3.2) 
where the Laplace transform of the memory kernel l|2.5|l reads as, 

_ ^ p jjggCjQ 

v ' p s Knn J_ 00 uj + Q- z v ; 
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Then from Eq. I|3.2(l we have, 

«.) - (3.4) 
— iz + V{z) 

and v(t) arises from the singularities of v(z) in the lower half-plane, Imz < 0. For instance, if the expression (|3.4|) has 
a unique simple pole located at zq = —W(zo), we get 

v(t) =v(0)e-* zot , (3.5) 

and the Markov approximation would be valid provided T)(zq) ~ D(0) (cf. Eq. I|2.7|l ). That is, taking the limit 
z — » i0 + in the Cauchy integral of Eq. 1|3.3[) we get^i 

P(0) = -±-P r + = d (3. 6 ) 

PsKIT J_ oc U + il p s K 

Therefore, from Eq. (|3.3(l we may realize that for such an approximation to be valid, it necessarily should be 
\ z o\ = M << CI, i- e -j CIDm(CI)/p s K << CI and ClD' M (Cl)/p s K << Cl. Now, according to the low-cyclotron frequency 
approximation (|3.1|l . the last condition will not be fulfilled for low enough frequencies, that is, the real part of zq will 
remain finite for Cl — > 0. This suggests the following approximation to find the poles of Eq. I|3.4|) : 

f>t \ -tVp n »(»-Rez ) p (°° uD AI (u)du (Cl - Rez ) 2 
iz = D(z ) ~ £>(Rez ) = — P / — — — 1 % fl - Rezo), (3.7) 

PsKTTll J _ co UJ + \l — HeZo p s K\l 



or, equivalently, 



I mZo = Re ^ 2 D M (Q-Rez ) (3.8) 



X 



(»-Rez ) p f uDM^du 2(Cl-Rez Q ) , , 
Rcz ° = FT - P / ~T^> — p o — / D M(u)du, 3.9) 

where the last equality arises from the approximation (|3.1|) . i.e., assuming £>m(0) ~ -DM(Oeff) << /q°° dui D m(w) / 'Q e ff> 
(f2 c ff = ^ — Rczo)- Then from Eq. I|3.9|) we get the solution 

/>oo 

fi off = Q - Rczo = Cl/{1 + [2/(p s KirCl)] / L> A /(w)dw}, (3.10) 

JO 

where, in fact, the effective frequency turns out to be bounded according to our previous discussion, < Cl c g < Cl. 
Note also that Cl e g - * corresponds to the limit of a vanishing cyclotron frequency, as expected. Finally, the friction 
coefficients D = — (p s K/f2)Irnzo and D 1 = (p s n/Cl)RezQ reads as, 

d = (ci cS /ci) 2 d m (ci cS ) (3.ii) 

D' = Ps k{1 - Cl eS /Ci), (3.12) 

which generalize the previous Markovian expressions H2.11|l and l|3.1|l . Here it is expedient to recall that (|3.11|) and 
(|3.12ll were extracted under the approximations of Eqs. i|3.1[) and (|3.7() . both being equivalent to D << D'. If, in 
addition, we have D' << p s K, then Cl e g ~ Cl and Eqs. I|3.11|l and (|3.12|) tend to the Markovian expressions. In other 
words, the frequency ratio Cl c g/Cl can be thought of as a measure of the proximity to the Markovian limit. Note that 
the limit of a vanishing cyclotron frequency is a strongly non-Markovian one, with D ~ 0(Cl 2 ) and p s K — D' ~ 0(0). 
Such a behavior of the transverse coefficient corresponds to the lower limit of Cl e g ~ 0(C1 2 ) as given by Eq. i|3.10[l . 
As regards the longitudinal coefficient, since the effective frequency was absent from our previous analysis in /, 
the limiting values for elastic scattering (CI — > 0) there reported are now drastically changed to vanishing values. 
However, we shall next see that there is a range of cyclotron frequency values that keep the previous agreement with 
the experimental values up to temperatures near 1 K. 



IV. ANALYSIS OF RESULTS FOR A PURE 4 He SYSTEM 
A. Longitudinal friction coefficient and the cyclotron frequency range 

In Table I we may compare values of the longitudinal friction coefficient computed from Eq. I|3.11|l for 17=0.01—0.03 
ps _1 , with the corresponding values arising from the Iordanskii formula 16 (phonon range), plus the R-R formula™^ 
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TABLE I: Longitudinal friction coefficient [10 -6 g cm -1 s -1 ] versus temperature for a pure 4 He system. The values in the 
third and fifth columns were calculated from Eq. I|3.1ip and have to be compared with the corresponding values in the second 
column, which arise from Refs. HH Il6l andl32t The values in the third and fourth columns were calculated for Q = 0.01 ps _1 , 
while the values in the fifth and sixth ones correspond to Q — 0.03 ps _1 . Powers of 10 are enclosed in brackets. 



1 [KJ 








M).03 


1 'Off / 1 ' 


0.1 


2.61[-8] 


2.63[-8] 


1.000 


2.97[-8] 


1.000 


0.2 


8.34[-7] 


8.15[-7] 


0.999 


8.45 [-7] 


1.000 


0.3 


6.33 [-6] 


5.98 [-6] 


0.996 


6.11[-6] 


0.999 


0.4 


2.84[-5] 


2.59[-5] 


0.989 


2.63 [-5] 


0.996 


0.5 


2.06[-4] 


2.07[-4] 


0.979 


1.95 [-4] 


0.993 


0.6 


2.43[-3] 


2.53[-3] 


0.962 


2.39 [-3] 


0.987 


0.7 


1.79 [-2] 


1.80[-2] 


0.938 


1.77[-2] 


0.979 


0.8 


8.27[-2] 


7.67[-2] 


0.900 


8.07[-2] 


0.964 


0.9 


0.274 


0.221 


0.840 


0.256 


0.940 


1.0 


0.714 


0.460 


0.749 


0.618 


0.900 



(roton range). Such calculations were performed taking into account only the phonon-roton contribution in Eq. I|2.11|) . 
i.e., for a pure 4 He system. In the phonon dominated regime (T < 0.4 K), the lack of experimental data leads us 
to compare with the results stemming from the Iordanskii formula, i.e., those arising from the limit ft — > of Eq. 
(|2.11|) . 17 Actually, it will be enough to have Uil/ksT < 1 in order to keep Djvf(O) ~ Dm(0) (cf. the figures for 
T = 0.1 K in Table I), but too small values of ft would affect the Markovian approximation ft c ff ~ ft, leading to 
appreciable discrepancies between the Iordanskii result, Dm(0), and our expression [j3.11Jl . The phonon influence 
becomes negligible for T > 0.65 K, allowing the comparison of our results with those arising from the R-R formula^ 
which constitutes a good fit to experimentally derived values up to temperatures about 1.3 Taking into account 
an estimated uncertainty of order 10 %, we may see from Table I that our results lie within such error bounds for 
temperatures up to 0.8 K (0.9 K) for ft = 0.01 ps" 1 (0.03 ps' 1 ). Note that for T=l K, our result for ft = 0.01 
ps _1 (0.03 ps _1 ) falls 30 % (4 %) below the lower error boundary. To summarize, we have identified a narrow 
range of cyclotron frequency values (0.01 ps" 1 < ft <0.03 ps _1 ) yielding longitudinal friction figures which are in 
agreement with the Iordanskii formula and experimental data. However, it is important to note that according to 
recent theories which show that the cyclotron frequency scales down logarithmically with the vortex size, such a 
range of ft would be consistent with a relatively wide range of macroscopic sizes of the system. 

B. Effective frequency, memory effects and transverse friction coefficient 

Another feature of Table I shows that the frequency ratio ft e ff/ft decreases with increasing temperature, which 
reflects a corresponding increase of D' (Eq. l|3.12[l L Particularly, at the lowest temperatures, the Markov approxima- 
tion ft e gf/ft ~ 1 shows to be excellent, becoming gradually less adequate as the temperature increases. The highest 
temperature range (0.9—1.0 K) displays the largest differences with the Markov approximation, as well as the highest 
discrepancy with the experimental results. In other words, the longitudinal friction phenomenon appears to be consis- 
tent with weakly, at most, non-Markovian processes (ft e ff/ft > 0.9). There is, however, another possible interpretation 
of such a discrepancy with the experimental results for ft e g/ft < 0.9, which is related to an eventual failure of the 
weak-coupling approximation. In fact, according to Eqs. I|3.10|l and (|2.11l) . we may see that for fixed ft, the parameter 
ft c ff/ft will behave as a decreasing function of the coupling strengths A and T, such that ft e s/ft — > 1 for a vanishing 
coupling (_Dj\/(cj) — ► 0), while ft c ff/ft — > for an infinite coupling (Dj{(aj) — > oo). Then, only a higher portion of the 
interval < ft e g/ft < 1 should be expected to be consistent with a weak-coupling approximation. In conclusion, the 
above discrepancy with the experimental results for ft e g/ft < 0.9 may also be regarded as an indication of a possible 
failure of the weak-coupling approximation. 

The transverse friction coefficient, in contrast to the longitudinal one, possesses a strong dependence on ft. In fact, 
Eq. (|3.1I) shows that in the Markovian limit, the bounds 0.01 ps _1 < ft <0.03 ps _1 lead to a factor 3 of spreading 
[D' M (0.0l ps -1 )=3-D^ f (0.03 ps -1 )], while the non-Markovian figures of Table I can reduce such a factor somewhat 
(>2.5). Another important difference between both friction coefficients stems from the degree of dependence on 
the quasiparticle dispersion relation cutoff featuresi^^i On the one hand, the longitudinal coefficient, which depends 
mainly on Dm(^hs)i turns out to be almost independent of such a cutoff, since ft ff and ft are two orders of magnitude 
lower than the roton frequency. The transverse coefficient, on the other hand, being mainly dependent on the integral 
J °° DM(<jj)dui, has therefore an important dependence on the quasiparticle cutoff through the corresponding depen- 
dence of the scattering amplitudes A and T, which is mostly uncertain. In summary, due to the above uncertainties in 
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the calculation of the transverse coefficient, only its order of magnitude should be reliable which, nevertheless, turns 
out to be quite useful to ensure that the condition D << D' is fulfilled. Recall that such a condition was assumed in 
the derivation of Eqs. Ij^.lljl and (|3.12(1 . and in fact, taking into account that p s n ~ 145 x 10~ 6 g cm~ 1 s _1 , all the 
values of Table I can be shown to be consistent with D « D' . It is worthwhile recalling also the lack of experimental 
results on D' for temperatures below 1.3 K. Taking into account only the vortex drag due to the scattering of rotons, 
the transverse coefficient can be written in terms of a transverse scattering length a±, viz. D' = p n VG<^±, where p n 
denotes the normal fluid density and vq the average group velocity of thermal rotons^*£ However, only speculative 
assumptions about the form of a± for temperatures below 1.3 K were reported^ In addition, it has been argued that 
the so-called Iordanskii forced gives rise to an additional transverse coefficient to be substracted from D', yielding a 
total transverse coefficientiii given by D t = D 1 — p n n. However, the sign, amplitude, and existence of this Iordanskii 
force are still subject to debatei^SSi^L^ Recently, Fortin^! has applied the formalism of Thouless, Ao, and NiuSi to 
compute the transverse and longitudinal coefficients due to the scattering of noninteracting phonons in two dimen- 
sions. He finds a transverse coefficient which turns out to be of opposite sign to ours and to that of Refs. and 01 
which is interpreted in terms of a negative vortex mass due to phonons. Such a discrepancy in the sign stems from 
the fact that, according to his equations, the transverse and longitudinal coefficients would be related, as functions of 
the cyclotron frequency, by Kramers-Kronig relations, while in our case such relations are connecting instead the real 
and imaginary parts of the Markovian frequency v (cf. Eq. H2.10|l '). 



C. Phonon emission and memory effects at zero temperature 



At this point, as a useful complement to our study, it will be instructive to discuss in some detail the memory 
effects related to phonon emission at zero temperature. We will base our analysis on the theory developed by Arovas 
and Freire^ for vortex dynamics in superfluid films. In fact, suppose that the vortex is set in motion at positive times 
by the action of a homogeneous time dependent superfluid flow, i.e., it is assumed that both, the superfluid velocity 
v s and the vortex velocity v, are zero for negative times. Then, the vortex equation of motion can be written^ 

t 

M(r)V(t - r)dr = i PsK [V(t) - V s (t)}. (4.1) 

o 

The right-hand side of this equation corresponds to the usual Magnus force p s nz x (v — v s ) (cf. Eq. expressed in 

complex notation (V = v x + iv y ), while the left-hand side will differ from the familiar Newtonian product of mass times 
acceleration, unless the memory or causal 4 kernel M(t) has a negligible lifetime. The memory, which actually plays 
an important role in this case, stems from the vortex coupling to the low lying excitations of the superfluid (phonons), 
in close analogy to the retardation and radiation effects stemming from electron-photon coupling in electrodynamics 4 . 
Now, we focus upon the long time limit of Eq. (|4.1|) . That is, for t >>lifetime of M(r), the left-hand side could 
be approximated as MV(t), where the effective vortex mass M is given by the Fourier (Laplace) transform of the 
memory kernel M(t) at zero frequency, M = drM(r) = M' + iM" , the imaginary part M" being related to the 
dissipation of vortex energy in the form of phonon emission. Then the solution of Eq. I|4.1|l in the long time limit 
could be easily obtained for constant V s : 

V(t) = V s {l-exp[i{n + iis r )t}}, (4.2) 
where the cyclotron Q and radiation damping v r frequencies respectively reads as, 

n = PsKM ' (43) 

M' 2 + M" 2 { ' 

- PbkM » 

r ~ A/' 2 + Af" 2 ' 1 ' 

M', however, is shown to diverge for an unbounded two-dimensional system 4 {M'{uj — > 0) ~ —In a;), while for a 
finite macroscopic system one should expect a corresponding finite value of M'(» \M"\), leading to the familiar 
expression <|1.3fl for the cyclotron frequency in l|4.3|l . with M' — m v the vortex mass per unit length. As regards the 
radiation damping frequency, from Eq. (14) of Ref. |j we have \M"(oj — ► 0)| = n 2 p s /{%c 2 a ), (c s = sound velocity) and 
then, v r — k(Q,/c s ) 2 /8. This result derives as well from the expression for the mean power radiated by unit length of 
a vortex performing cyclotron motion (see Eq. (2.11) of Ref. 0). Therefore, the radiation damping should be weak 
v r ~ 10~ 3 r2 for our range of cyclotron frequency values (f2 ~ 0.01 ps~ 4 ), whereas it would be relatively strong, v r ~ Q, 
for cyclotron frequencies arising from a hydrodynamical model for the vortex mass^ fl ~ 3 ps _1 . It is interesting 
to notice that a vanishing damping frequency for an infinite system in (14 .4|) does not preclude the approach of the 
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TABLE II: Longitudinal friction coefficient [10 _6 g cm -1 s _1 ] versus temperature and 3 He concentration for mixtures. The 
third column corresponds to experimental results, while the forth one corresponds to values arising from Iordanskii formula 1 — 
(phonon contribution) plus R-R formulas 14 (roton+ 3 He contribution). The remaining notation is the same as in Table I. 



T[K] 


C 


^exp 


^Rcfs 


Dq.01 


n c ff /n 






A). 03 




0.28 


2.84[-5] 


4.69[-3] 


4.69[-3] 


4.69[-3] 


0.995 


4.69[-3] 


0.997 


4.64[-3] 


0.998 


0.28 


7.55[-6] 


1.35[-3] 


1.25 [-3] 


1.25 [-3] 


0.996 


1.25[-3] 


0.998 


1.24[-3] 


0.999 


0.61 


7.55[-6] 


4.56[-3] 


5.01 [-3] 


4.87[-3] 


0.960 


4.90[-3] 


0.980 


4.80[-3] 


0.986 


0.615 


1.4[-7] 


3.34[-3] 


3.43[-3] 


3.56[-3] 


0.959 


3.52[-3] 


0.979 


3.39 [-3] 


0.986 


0.643 


1.4[-7] 


6.04[-3] 


6.17[-3] 


6.36 [-3] 


0.953 


6.32[-3] 


0.976 


6.11[-3] 


0.984 


0.67 


1.4[-7] 


1.01 [-2] 


1.04[-2] 


1.07[-2] 


0.947 


1.07[-2] 


0.973 


1.04[-2] 


0.982 



vortex velocity to the superfiuid velocity at long times. Actually, it simply means that such an approach will be slower 
than the exponential one of Eq. (|4.2[1 . To see this, let us integrate by parts the left-hand side of Eq. (|4.1|l getting 
M(0)V(t) + f M(r)V{t — r)dr. Then, approximating in the long time limit the last integral as f Q M(t)V (t)dr , the 
left-hand side of Eq. (|4.1|l turns out to be simply M(i)V(t), from which we get 



V(t) = V s 



1 - iM(t)/p s K 



1 + {M{t)/ Ps Kf 



(4.5) 



wheroi M(t)/p 8 K ~ £/(2c s t), £ — n/(2irc s ) being the coherence length. Thus, the above expression shows that in the 
case of an infinite system, the approach of the vortex velocity to the superfiuid velocity turns out to be, in contrast 
to l|4.2[) . a slow nonexponential one. 



V. ANALYSIS OF RESULTS FOR DILUTE SOLUTIONS OF 3 He IN 4 He 

In case of a 3 He- 4 He mixture we have to take into account both terms in the expression (|2.11|) for Dm(u>). The 
calculation of the fermion term D 3 (uj) (cf. Eq. Q2.12)) ) turns out to be similar to that leading to D 3 (0) in Eq. i|2.13[l . 
namely 

* A2 poo 

D ^ = ^^ dk\T(k,q)\ 2 (f k -f q )k 2 (<f+k 2 /3), (5.1) 
it Ti ui Jo 

where the value of the momentum q arises from the argument of the second Dirac delta in Eq. H2.11[) . i- e -> q 2 = 
k 2 + 2m*ui/h. Thus, the cutoff of the Landau- Pomeranchuk dispersion relation imposes the same cutoff (~ 1 ps _1 , 
see Ref. l33|) to the frequency lo in Eq. (|5.1|l . Recall that the cutoff uncertainties will be reflected in the evaluation 
of the transverse friction coefficient, as was already mentioned in Sec. IIVBI To compute Eq. I|5.1|) we utilized the 
following simple generalization of the expression (|2.14|l for k ^ q: 

128 (m*Ay 

In Table II we may compare some experimental results for the longitudinal friction coefficient (third column) , with 
the corresponding results arising from our approach, along with the values computed from the Iordanskii 16 and R-R 14 
formulas. Given the low 3 He concentrations of Table II, the Fermi temperatures turn out to be at most two orders of 
magnitude below the experimental ones, and so the Fermi occupation numbers in Eqs. I|2.13|l and (|5.1I) can be very well 
approximated by the Maxwell-Boltzmann statistics. At T=0.28 K the phonon contribution to the friction, stemming 
from the Iordanskii formula, turns out to be negligible in comparison to the 3 He term given by Eqs. (|2.13(l and (|2.14(l 
(cf. the values of -Dpiofs in Tables I and II). Then, the experimental value D oxp for T=0.28 K and C=2.84xl0~ 5 in 
Table II (being C = na/foa + 714), rii= number density of ^e atoms), which was measured within an error less than 
1%, was utilized by R-R to calculate the effective cross section cto = 18.3 A in Eq. H2.14|l . assuming for the effective 
mass the value m* = 2.5m3 (m^= actual mass of a 3 He atom). A similar procedure was followed in our case, since 
the value of ao in i|5.2f) was set to get agreement with the experimental value Z? exp =4.69xl0 -3 for fl = 0.02 ps _1 , 
i.e., the center of the cyclotron frequency range discussed in Section llV Al Then, equating (|3.11|l to 4.69xl0~ 3 , we 
extracted the value Co = 18.54 A, which turns out to be slightly greater than the R-R result. Although the phonon 
contribution to Dm in Eq. (|3.11|) is in fact negligible for T=0.28 K and C=2.84xl0~ 5 , this is not the case for Sl c ff/Sl 
(Eq. (|3.10|l ). since the phonon contribution to the integral L°° DM{u)duj turns out to be greater than the 3 He one. 
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Nevertheless, as seen from Table II, the factor (f2 e ff/f2) 2 in (|3.11() remains close to unity and then, in practice, almost 
all the longitudinal friction should be ascribed to 3 He scattering. 

In a second experiment, to prove the proportionality of the friction coefficient to the 3 He number density in the dilute 
limit, R-R performed a measure for the same temperature T=0.28 K and a smaller concentration of C = 7.55 x 10 -6 . 
They obtained the value D cxp — (1.35 ± 0.06) x 10 -3 , which turns out to be in agreement with their theoretical 
calculation arising from Q2.13(l . within the limits of estimated error. As regards our calculation from Eq. H3.11J) . it 
yields practically the same figures as the R-R formula (Table II). 

The diluted sample was also used to verify the additivity of 3 He and roton scattering, by performing an experiment 
at the relatively high temperature of 0.61 K. The directly measured value D exp — 4.56 x 10~ 3 was then contrasted 
with that arising from the addition of R-R formulas for 3 He and roton scattering contributions, namely 4.79xl0 -3 . 
The phonon contribution, on the other hand, was ignored, presumably because of some discrepancies arisen from the 
Pitaevskii'siii calculation of the phonon-scattering cross section. Actually, taking into account such a contribution, the 
value of the friction coefficient should have been increased to 5.15xl0~ 3 . The Pitaevskii's result was later modified 
by Iordanski: 16 > in that the coefficient of proportionality to T 5 of the friction coefficient due to phonons was shown to 
be smaller by a factor ~ 0.62. That is, taking into account the phonon contribution stemming from the Iordanskii 
formula, the corrected value 5.01 xl0~ 3 (Table II) is in fact closer to the experimental one. Finally, we compare with 
our results computed from Eq. (|3.11() . From Table II we see that such results are closer to the observed one than the 
previous estimation of Iordanskii+R-R, and in this better agreement it is important to remark the role played by the 
memory effect, which is embodied in the factor (n c ff/tt) 2 < 1 in Eq. 13.11fl . 

Next we analyze a set of measures performed for ordinary helium (C = 1.4 x 10 -7 ) at temperatures of 0.615, 0.643 
and 0.67 K. The second measure (T = 0.643 K) was reported in Ref. H], while the remaining two ones are included 
in Ref. ^3- Under such conditions, the incidence of 3 He scattering is almost negligible in both, the Iordanskii+R-R 
results and our figures computed from (|3.11|) . From Table II we see that, analogously to the above results for T=0.61 
K, theoretical calculations again overestimate the experimental data, and the best agreement is also obtained for our 
result at Q — 0.03 ps _1 . 

Finally, from Table II we notice that the memory effect remains small (f2 c ff/f2 > 0.95), showing the same increase 
with temperature as in Table I. On the other hand, the dependence of Q c g/Q on concentration is not clear from Table 
II except for T=0.28 K, where we find a slight reduction for a higher concentration. It is not difficult, however, to 
generalize such a result taking into account that in the dilute limit, the Maxwell-Boltzmann approximation to fk and 
f q in Eq. (|5.1|l yield a D^(lo) proportional to the 3 He number density, which in turn implies a growing of DM{u)du> 
with concentration and, accordingly, a decreasing behavior for il c ff/S7 in Eq. I|3.10|l . In other words, f2 e ff/f2 is shown 
to be a decreasing function of the number of quasiparticles, i.e., phonons, rotons and 3 He impurities, interacting with 
the vortex. 

We hope, in concluding this report, that the present theoretical results will stimulate an experimental investigation 
of degenerate- 3 He richer vortex friction regimes. 



Our starting point is the Hamilton equation of motion for the creation operator of right circular quanta, whose 
time derivative turns out to be proportional to the vortex velocity: 



APPENDIX: DERIVATION OF THE VORTEX EQUATION OF MOTION 



a) = \J npsL/m^R, 



(A.l) 



where L denotes the vortex line length. Thus we have, 




(A.2) 



where O denotes the quasiparticle operator: 



O = 



^2 5 Kq z [{ k v - %) +i{(lx - fc x )][A(/c,o)6j c 6 q + T(k,q)cl i 




(A.3) 



Next we consider the Hamilton equation for d^: 




^ fa,qA(K ~<ly) + i (Qx - k x )][A(k,q)(u k 
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+ T^(k,q)(ek -e 9 )4<r c q^] ~ 72 — ~r, Yl Yl ^k'Jk'.qA^x - K)i k ' v ~ k y) 

+(k' y - q y )(k' x - k x )][A(k, k')A(k', q)bt 6 q + T(k, k r )T(k', q)cl a . (A.4) 

As in /, we shall make use of a weak-coupling approximation, which consists in retaining only second order terms in 
the scattering amplitudes A and T. This approximation is discussed in Sec. IIVBI Then to approximate the above 
expression, we note that to the zeroth order in such parameters we have 



a 



-ia</Sl (A.5) 

&£(<)& q (t) = e*^*-"*'* 6^(0) 6,(0) (A.6) 

4,a(*)Cq.«r(t) = ^-^^4^(0)^(0), (A.7) 

while the first order arises from the Hamilton equations for 6 k 6 q and 4 a c qi(T . We have, for instance, 

= i(u k - u q )bl 6 q + * {^2 S Kkz A(k', k)[[(k y - k') + i(k x - k' x )]a ] 

hyJ-Kp s L/m 4 k , 

+[(k' y - k y )+i(k x - k' x )]a] b\. t 6 q - ^ S q ^ A(q, q')[[{q' y - q y ) + i(q' x - q x )W 

q' 

+[(qv - iy) + »(«£, - ?*)]«*] 4 M> ( A - 8 ) 

and we may find a formal solution to this equation in b^(t) 6 q (i) by noting that the only dependence on 4 6 q on the 
right-hand side comes from the first term. Thus we have to the first order in A: 



4(i)6,(i) = e^-^)*6] £ (0)6 q (0)+ / dre*^*— ====(£ A(fc', k)[[{k v - k' y ) 

Jo hy / irp s L/m 4 k , 

+i{k x - k' x )]h\t - r)/n + [(&„ - fc;) + - fc x )]d(t - r)/fi] e^'-^-^bliO) 6,(0) 

_ <V„5. A (<?'9')[[(^ - «v) + - Qx)]a j (t - r)/fi + - g v ) 
q' 

+<(& - <£)]d(t - r)/fi] e i ^-^')(*-^6 k (0) 6,,(0)}, (A.9) 

and analogously we may find a similar expression for 4 a (t) c q . cr (t). Finally, replacing Eq. (|A.9Jl (and the corre- 
sponding expression for c t a (t) c, )CT (f)) in Eq. I|A.4(> and taking mean values according to (6 k (0) 6,(0)) = S^qUk and 
(4 <r(°) c q^( )) = ^kq/fc, we obtain Eq. 
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